Transport and localization in periodic and disordered graphene superlattices. 
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We study charge transport in one- dimensional graphene superlattices created by applying layered 
periodic and disordered potentials. It is shown that the transport and spectral properties of such 
structures are strongly anisotropic. In the direction perpendicular to the layers, the eigenstates in 
a disordered sample are delocalized for all energies and provide a minimum non-zero conductivity, 
which cannot be destroyed by disorder, no matter how strong this is. However, along with extended 
states, there exist discrete sets of angles and energies with exponentially localized eigenfunctions 
(disorder-induced resonances). Owing to these features, such samples could be used as building 
blocks in tunable electronic circuits. It is shown that, depending on the type of the unperturbed 
system, the disorder could either suppress or enhance the transmission. Remarkable properties of the 
transmission have been found in graphene systems built of alternating p-n and n-p junctions. The 
mean transmission coefficient has anomalously narrow angular spectrum, practically independent of 
the amplitude of the fluctuations of the potential. To better understand the physical implications 
of the results presented here, most of these have been compared with the results for analogous 
electromagnetic wave systems. Along with similarities, a number of quite surprising differences 
have been found. 



PACS numbers: 71.10.Fd, 73.20.Fz, 73.23.-b 
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INTRODUCTION. 



The exploration of graphene is nowadays one of the 
most animated areas of research in condensed matter 
physics (see, e.g., Refs.^^^^^). Its unique properties not 
only arouse pure scientific curiosity but also suggest pos- 
sible practical applications. More in-depth studies of 
graphene continuously bring about more counterintuitive 
discoveries. Examples are plentiful. Suffice to men- 
tion a novel integer quantum Hall effect^^^, total trans- 
parency of any potential barrier for normally-incident 
electrons/holes^ (in analogy with the Klein paradox^), 
and the recently predicted focusing of electron flows by a 
rectangular potential barrier^ (an analog of the Veselago 
lens^ii^). Of even greater surprise are the properties of 
disordered graphene systems -^-^^ ^-^^i-^^i-^^i-^^i-^^i-^^i-^^ . The lat- 
est results, both theoretical and experimental, led to the 
amazing conclusion that there is no localization in dis- 
ordered graphene, even in the one-dimensional situation, 
i.e., when the random potential depends only on one co- 
ordinate. In this paper, we show that this conclusion, if 
taken unreserved, could be misleading. We demonstrate 
here that a well-pronounced localization can take place in 
graphene, i.e., there could exist a (quasi)-discrete spec- 
trum with exponentially localized eigenfunctions^^. This 
localization can occur even though disorder can never 
make a graphene sample a complete insulator, and there 
is always a minimal residual conductivity (an indication 
of delocalization). In this paper, the charge transport in 
periodically and randomly layered graphene structures is 



studied and analogies with the propagation of light in 
layered dielectrics are discussed. 



II. BASIC EQUATIONS 
A. Charge transport in graphene 

A graphene layer consists of two triangular sublattices 
(A and B). The low-energy band is gapless and electronic 
states are found near two (electron and hole) cones. The 
behavior of charge carriers (electrons and holes) near the 
Dirac point is governed by the 2D Dirac equatio n^^i^^ : 

vf(^-p)^ = (^- V)^, (1) 

where ^ is a two-component spinor (^a, ^b)"^, the com- 
ponents of a pseudospin matrix a are given by Pauli's 
matrices, p is the momentum operator, vf is the Fermi 
velocity, V{x^ y) is the potential, and E is the state en- 
ergy. When the potential depends on one coordinate, 
V = V{x)^ the wave function ^f{x^y) can be written as 
^{x^y) = e*^^^V^(x), and Eq. ^ can be presented in the 
dimensionless form 
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Here ^ = x/d^ d is the characteristic spatial scale of the 
potential variations, e = Ed/Tivp^ u = Vd/Tivp^ and 
P — kyd. 
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In what follows, we consider potentials ix(^) comprised 
of periodic or random chains of rectangular barriers de- 
picted in Fig.[T]In a j-th layer, the solution of Eq. ([2|) has 

























1 

1 















FIG. 1: Schematic diagram of the potential u{^). 



The spinor components at the left and right boundaries 
of the j-th layer are connected by the diagonal matrix 
Sj = diag(e*"j , e~*"j ), where aj = hzjdj is the phase ac- 
cumulated by the wave propagating through the layer of 
the thickness Sj. Thus, the matrix Sj-^iAdjj-^i transports 
the spinor components from the left side of the interface 
between the j-th and {j + l)-th layers to the left side of 
the next interface between the {j + l)-th and {j + 2)-th 
layers (see Fig. [2j). Obviously, the total transfer matrix 
of a layered sample consisting of N layers is given by the 
product of SM matrices: 

N 

f=n4+i^w+i (8) 

i=o 



the for m: ,p^A,Bh = ^ (a,b)/"''^ + ^(a,B), 
Kj = amplitudes of the 

rightward (+) and leftward (— ) propagating spinor com- 
ponents. At the layer interfaces the amplitudes V^^^^^^ . 
and V^l^^^) .^^ are connected by the equation 



^, where 
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(3) 



which follows from the continuity of the spinor compo- 
nents. Since from Eq. ([2]) the amplitudes '^[^^^^ . are 
connected. 



FIG. 2: (color online) Propagation through a layered struc- 
ture. The matrix Sj propagates through j—ih layer (green 
arrow) and the matrix Mjj-^i connects the spinor amplitudes 
across the interface between j-th and {j + l)-th layers. 



e — Uj 



(4) 



B. Light transport in dielectrics 



= o Sg^(g-^j) = ±e ' •sgn(£-ix^-). (5) 

we will only consider the amplitudes ?/^^^ and omit 
the subscript "A". If \m.Kj = 0, 6j = arctan(/3/A>:j) is 
the angle between the wave vector k of the (+) wave 
and the normal to the interface between the j-th and 
{j -\- l)-th layers (i.e., the angle of propagation in the 
j-th layer). Using Eqs. (|3]) and (|4]) one can calculate 

the matrix Mjj-\-i that connects the amplitudes ?/;-^^ 



and V^j+i on two sides of the interface, ^V^j+i, V^j+i 
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2 cos^j+1 
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where 



= e-''^+^ ± e±'^^- • sgn[(e - u,){e - u,+^)], (7) 
and T denotes transposed vector. 



Analogous products of matrices have been well stud- 
ied in the context of transport of electromagnetic waves 
in layered media (see, for example, and references 
therein). To better understand the physics of charge 
transport in graphene subject to a coordinate-dependent 
potential, in what follows, we contrast the results for 
graphene with those for the propagation of light in lay- 
ered dielectric media (for more analogies between quan- 
tum and optical systems, see, e.g. Ref.^^^^^). Additional 
analogies, not discussed here, also exist with the trans- 
port and localization of phonons in different kinds of pe- 
riodic and random one-dimensional structures^^^^^^^^. 

In the latter case, the matrices Sj are the same as in 
Eq. (|8]), and the transfer matrix, A^j^^+i, 
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(9) 



that describes the transformation of the amplitudes of 
the electromagnetic waves at the interface between jth 
and (j + l)th layers, has the form Eq. (|6]) with gj^\_i 
being replaced by 



cos6>j+i ± cosOj • sgn(njnj+i)- 



(10) 
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for 5-polarized waves and 



cos ± COS • sgn(njnj+i) (11) 



for p-polarized waves. Here, is the angle of the propa- 
gation, Zj = ^J jiij I Sj is the impedance of jth layer, and 
rij = -^yfejjij is its refractive index. The signs ± cor- 
respond, respectively, to dielectrics with positive (right- 
handed, R) and negative (left-handed, L) refractive in- 
dices. 

It is easy to see that the parameter {e — u) plays, in 
graphene, the same role as the refractive index n in a 
dielectric medium. It is due to this similarity that a p- 
n junction (interface between regions where the values 
{e — u) have opposite signs) focuses charge carriers in 
graphene, like an R-L interface focuses electromagnetic 
waves ^. 

Note that in Eq. (O (for graphene) there is no fac- 
tor Zj+i/Zj, which determines the reflection coefficients 
at the boundary between two dielectrics^^. This means 
that the charge transport in graphene is similar to the 
propagation of light in a stack of dielectric layers with 
equal impedances. In particular, both p-n and p-p 
junctions are transparent for normally incident charged 
particles^'^. This property is readily seen from the anal- 
ysis of the matrix M: at /3 = it is a Pauli matrix 
M = (Tx for p-n junctions and unit matrix for p-p junc- 
tions. Therefore, a (+)-wave is totally transformed into 
a (— )-wave at a p-n junction, and remains a (+)-wave 
at a p-p junction. Another important difference between 
the transfer matrices M (graphene) and Ai (electromag- 
netic waves) is that M is, generally speaking, a complex- 
valued matrix, while the is always real. As it is shown 
bellow, this distinction brings about rather peculiar dis- 
similarities between the conductivity of graphene and the 
transparency of dielectrics. 



III. TRANSPORT IN PERIODIC STRUCTURES 

Among the vast amount of publications on graphene, 
a significant and ever increasing part belongs to pa- 
pers devoted to the charge transport in graphene su- 
perlattices formed by a periodic external potential (see. 



e.g.f 



:). This is not only due to its theoret- 
ical interest but also because of the possibility of ex- 
perimental realization and potential applications^. In 
Ref.^, for example, it was suggested that, by virtue 
of the high anisotropy of the propagation of carriers 
through graphene subjected to a Kronig-Penney type 
periodic potential, such a structure could be used for 
building graphene electronic circuits from appropriately 
engineered periodic surface patterns. 

Here we consider a layer of graphene under a periodic 
alternating potential u = ^Uq and assume that 6 = 0. 
Two layers, ui = Uq and U2 — —Uq with thicknesses 
6i and ^2, respectively, constitute the superlatice period 



with transfer matrix T = /S'2M2,i'S'iMi,2- Its eigenvalues 
X{/3) indicate whether this periodic structure is transpar- 
ent or not. Namely, the structure is transparent when 
\X{P)\ = 1 and is opaque when |A(/3)| 7^ 1. Note that 
analogous (ni = —712) L-R periodic dielectric structures 
are transparent at all angles of incidence when Zi = Z2. 
In contrast, a periodic array of p-n junctions in graphene 
has a rather nontrivial angular dependence of the trans- 
mission coefficient T{P). This distinction between pe- 
riodic graphene and dielectric lattices follows from the 
difference in the corresponding transfer matrices: M is 
complex- valued while Ai is real. 

If £ = and Si = 62 = S (symmetric graphene system) 
the equation for the eigenvalues of the matrix T has the 
form 



-2A 



1 — sin^ cos 2a 
cos^ 6 



0, 



(12) 



where a = hz{0)6. It is easy to see that |A| = 1 at normal 
incidence {0 = 0) and at a discrete set of angles, Om, 
given by 



a{0m)=Tn7r^ m = ±l,±2,... 



(13) 



Although the eigenvalues A(/?) are well defined for an 
infinite system, they are also quite meaningful for a suffi- 
ciently long finite periodic sample: at Om there are max- 
ima of the transmission coefficient T{[3). The transmis- 
sion coefficient T{[3) of the symmetric graphene structure 
is presented in Fig. [3^. a similar transmission spectrum 
T(/3) exists at 6> 7^ in L-R periodic structures (Fig. [3)3) 
made of layers with equal absolute values of the refrac- 
tive indexes, ni = — n2, and different impedances Zi ^ 
^^35,36 rpj^^ transmission coefficient T(/3) for a symmet- 
ric periodic R-R system is shown in Fig. [3b. 

The transmission spectrum T{(3) in Fig. [3^, which con- 
sists of a discrete set of incidence angles, is the result of 
the degeneracy caused by the high symmetry of the struc- 
ture {ui — —U2^ e = 0, and Si = S2). Any symmetry- 
breaking splits the degeneracy and the spectrum takes 
the form usual for ordinary periodic structures: a set of 
conducting zones of non-zero width separated by band 
gaps^^. In Fig. 2^, the zone structure of the transmis- 
sion spectrum is shown. It is important to note that 
instead of the wave number ky and the energy, typically 
used in zone diagrams, the variables in Fig. [4^ are the 
asymmetry parameter A = {Si — 62) /2 and P = kd sin 0, 
respectively. Note that even in non-symmetric structures 
there are some values of Si and S2 for which, along with 
the usual conducting zones and band gaps, there exists a 
discrete set of resonant /3s. Note also that, for a fixed A, 
the transmission zones as a function of P are very narrow, 
making the direction of the charge flux easily tunable by 
changing the applied voltage Uq. 

For comparison, the analogous spectra for L-R and R- 
R periodic structures with \ni\ = \n2\ and Zi 7^ Z2 are 
presented in Fig. [4)3, c. One can see there that the trans- 
mission coefficients of graphene and L-R structures are 
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FIG. 3: Transmission coefficient T{f3) for: (a) grapfiene sub- 
ject to a symmetric periodic potentiai u = itt/o; (b) a sym- 
metric periodic L-R dieiectric structure; (c) a symmetric pe- 
riodic R-R dielectric structure. Tlie parameter Uo plays in 
graphene the same role as the refractive index no in a dielec- 
tric medium. For (b) and (c), no = Uo and Z1/Z2 — 1.1. In all 
figures, L stands for left-handed media and R for right-handed 
media. Thus, the structures considered in (b) are periodic 
stacks LRLRLR. . . of dielectrics. The structures considered 
in (c) are periodic RRRRRR. . . stacks of dielectrics. 



similar, and both differ drastically from that of the R-R 
structure. 

A phenomenon similar to the total internal reflection of 
light can occur to charge in graphene at a non-symmetric 
p-n interface when \e — ui \ 7^ \e — U2\ (here ui and U2 are 
the potentials on either side of the interface). However, 
a periodic set of such junctions is transparent for some 
angles of incidence^^. This effect is similar to photon 
tunneling (frustrated total internal reflection) in stacks 
of dielectric layers^^. 

Let us now consider the transmission of charge through 
a single non-symmetric p-n junction, assuming that Ui = 
-U2 = Uo > e > 0. Then, if {e - U^f < /^^ < (e + 
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FIG. 4: (color online) Transmission coefficient T[l3^ A) as a 
function of the difference of the thicknesses of two layers, A = 
(^1 — ^2)72, and of the normalized parameter jS — kyd = 
kd sin for: (a) periodic graphene subject to an alternating 
periodic potential; (b) periodic L-R dielectric structure; (c) 
periodic R-R dielectric structure. In (a), the normalization 
parameter Uo is the dimensionless applied voltage; in (b) and 
(c), the normalization parameter no is the refraction index 
in a dielectric medium, and Z1/Z2 = 1.1. The grey area 
corresponds to a perfect transmission, T — 1. The white 
regions correspond to T = 0. 



/7o)^7 a total internal reflection occurs at the interface; 
however there are ranges of the angle of incidence where 
the system is transparent. An example of such an angular 
spectrum is shown in Fig. [5^. 

Physically, the tunneling conductance (proportional to 
the transmission coefficient) is due to the confined states 
in graphene quantum well^^^^^. Although the confined 
states in a single quantum well have a discrete spectrum 
P = Pni^)^ an infinite periodic chain of wells, interacting 
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FIG. 5: (color online) Transmission coefficient T{l3) for: (a) 
periodic graphene subject to an alternating potential; (b) pe- 
riodic L-R dielectric structure; (c) periodic R-R dielectric 
structure. The coordinates and the colors are the same as 
m Fig. H In the region above the straight line the condition 
of total internal reflection is satisfied. 



via their evanescent wave functions, forms transmission 
bands centered around jSn- 



Periodic L-R and R-R dielectric structures have prop- 
erties similar to the properties described in this subsec- 
tion for graphene. Figures [5)3, c show the transmission 
spectrum of L-R and R-R structures with a period com- 
posed of two blocks with equal thicknesses and different 
refractive indexes, |ni| = 1 + (5n and |n2| = 1 — 5n. 



IV. TRANSMISSION IN DISORDERED 
STRUCTURES 



Based on the results obtained for ideally periodic sys- 
tems, the authors of Ref. [sllssi suggested that graphene 
superlattices could be used as tunable elements in elec- 
tronic devices. Since parameters of such structures are 
extremely sensitive to the variations of the applied po- 
tential it is worthwhile to study the effect of disorder 
(random deviations of the potential from periodicity) on 
the propagation of charge in such configurations. More- 
over, this study is of interest by itself because strongly 
disordered (with no periodic component) potentials bring 
about further unexpected spectral and transport proper- 
ties of graphene samples, which make them potentially 
useful as an alternative to pure periodic systems. 

A surprising and counter-intuitive result is that a sam- 
ple of graphene subject to a random one-dimensional po- 
tential, u{x)^ is absolutely transparent to the charge flow 
perpendicular to the x-direction, no matter how long is 
the sample and how strong the disorder isi^. This means 
that in such samples there exists a minimal non-zero con- 
ductivity, which (together with symmetry and spectral 
flow arguments) led to the conclusion that there is no lo- 
calization in ID disordered graphene systems^^^^^. How- 
ever, this statement (being correct in some sense) should 
be perceived with a certain caution. Below, we show that 
although the wave functions of normally incident parti- 
cles are extended and belong to the continuous part of 
the spectrum, away from some vicinity of 6> = 0, 1-D ran- 
dom graphene systems manifest all features of disorder- 
induced strong localization. Indeed, there exist a discrete 
random set of angles (or a discrete random set of energies 
for each given angle) for which the corresponding wave 
functions are exponentially localized with a Lyapunov ex- 
ponent (inverse localization length) proportional to the 
strength of the disorder. 

Obviously, the behavior of a quantum mechanical par- 
ticle is determined by the type of potential and by the 
ratio between its values for u{x) and the energy e of 
the particle. Bellow, we study the charge transport in 
graphene subject to a random layered potential of the 
form Uj = UQ{j) + AiXj. Three particular cases are con- 
sidered: (i) all Uj < £, uo{j) is a periodic function; 

(ii) £ < uo{j) = const; and (iii) e = 0, uo{j) is a 
periodic set of numbers with alternating signs. In all 
cases, the Auj are independent random variables homo- 
geneously distributed in the interval [—Auq^Auq]. In 

(iii) , \Auj\ < \uo{j)\ and uo{j) = it/o, which represents 
an array of random p-n junctions, where electrons out- 
side a barrier transform into holes inside it, or vice versa. 
For the sake of simplicity, here we assume that the widths 
of the layers do not fluctuate, will be addressed latter. 
These three cases will be considered in the next three 
subsections. 
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A. Case (i): all Uj < s, and uo{j) periodic 

Figure Fig. [6] shows an example of the angular depen- 
dence of the transmission coefficient, T(/?), for a type (i) 
graphene sample that contains 40 layers of equal thick- 
ness ^0 = 1-0, ui — —7.0, U2 = —13.0, e = 0.0. One 
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FIG. 6: (color online) Transmission coefficient T(/3) for 
periodic (thin black line) and disordered (bold blue line) 
graphene. The range of the variation of the potentials, 
Auo O.lUo; s = 0. 

can see that a relatively weak disorder has drastically 
changed the transmission spectrum: all features of the 
spectrum of the underlying periodic structure has been 
washed out, and a rather dense (quasi-) discrete angu- 
lar spectrum has appeared with the corresponding wave 
functions localized at random points inside the sample 
(disorder-induced resonances). Figure [71 shows the spa- 
tial distribution of the square modulus of the amplitude 
of a resonant wave function (intensity distribution inside 
the sample). For a fixed £, T(/3), shown in Fig. [6l has 
the same form as T{£) for a fixed /3, shown in Fig. [51 
Both consist of randomly distributed resonances (one in 
the f3 domain and another versus energy) typical for ID 
Anderson localization of electrons and light. However, 
there is one fundamental difference from the usual An- 
derson localization: in the vicinity of normal incidence, 
the transmission spectrum of graphene is continuous with 
extended wave functions, and the transmission coefficient 
is finite (T = 1 at ^ = 0). It is this range of angles that 
provides the finite minimal conductivity, which is pro- 
portional to the integral of T{6) over all angles 0. 

The mean transmission coefficient, {T{/3)), for different 
strengths of disorder (different Ai^o) is plotted in Fig. [9l 
As expected, the increase of disorder reduces the trans- 
mission and narrows down the angular width of the trans- 
mission spectrum A/3. The zero (to within the resolution 
of the plots) values of {T{/3)) at each curve correspond 
to the angles^ which exceed the angle of total internal 




FIG. 7: Spatial distribution of the wave function localized 
inside the sample for P marked by a red arrow in Fig. [6] 
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FIG. 8: (color online) Transmission coefficient as a function 
of the energy £, T{s), for periodic (thin black line) and disor- 
dered (bold blue line) graphene. Ai^o = 0.1[/o, f3 = O.SUq. 



reflection. 



B. 



Case (ii): e < uo{j) = const 



In this case, the results are more intriguing^^-i^ (al- 
though encountered in usual electron and optical ran- 
dom systems^^). In this case, the transmission of the 
unperturbed system is exponentially small (tunneling) 
and gets enhanced by the fluctuation of the potential 
(Fig. [To]). This is quite natural because the transmission 
of each j-th segment is proportional to ex.p[—6o{e±Auj)] 
and (despite of the fact that (Auj) = 0) the mean value 
(exp[— (5o(£ ± Aiij)]) > exp(— (^0^)7 due to the asymme- 
try of the exponential function. Note that another type 



7 



Disordered graphene 




FIG. 9: Mean transmission coefficient {T((3)) for disordered 
grapfiene. Tfie curve (1) corresponds to Auq/Uq = 0.1, tfie 
curve (2) corresponds to Auq/Uq = 0.2, etc. Curves (8), (9), 
and (10), for wfiich Auq/Uq = 0.8, 0.9, and 1.0 are practicaiiy 
indiscernibie. 
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FIG. 10: Mean transmission coefficient {T(f3)) for disordered 
graphene for different vaiues of the parameter Auo/Uo, which 
increases from zero (narrowest curve) to 1.0 (widest curve). 
Here, Auq/Uq increases by 0.1 from ieft to right. 



of disorder, linked to graphene layer edges, leads to the 
same result: the disorder improves the transmissioni^, 
compared to the ordered graphene case. 



C. Case (iii): \Auj \ < \uo{j)\ and uo{j) = ±Uo 

The behavior of the charge carriers in the graphene 
system of type (iii) is most unusual. It is characteris- 
tic of two-dimensional Fermions and have no analogies 
in electron and light transport. Shown in Fig. [11] by the 
bold blue line is the transmission spectrum at £ = of a 



graphene sample containing 40 layers of equal thicknesses 
Sj = 5o and alternating random potential. One can see 
that, compared to the underlying periodic configuration 
(thin black line), the disorder: obliterates the transmis- 
sion peaks located near with m 7^ 1 [see Eq. (p!3|) ]: 
makes much wider the transparency zone near /3 = 0; 
and gives rise to a new narrow peak in the transmission 
coefficient, associated with wave localization in the ran- 
dom potential. In contrast to the peaks in the periodic 
structure, the wave function of this disorder-induced res- 
onance is exponentially localized. 



T{/]) 




Disordered graphene 



0.4 0.6 



0.8 



1.0 



FIG. 11: (color online) Transmission coefficient T{f3) for peri- 
odic (thin black line) and random (bold blue line) symmetric 
potential with Auj G (—0.1,0.1). 
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FIG. 12: (color online) Mean transmission coefficient {T(f3)) 
for disordered graphene. Arrows with numbers mark the 
strength of the disorder, Auq/Uq. 

For this case (iii), the average transmission coefficient 
as a function of /3 is presented in Fig. [121 In contrast to 
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the case (i), the transmission in (iii) is extremely sensi- 
tive to fluctuations of the applied potential: in Fig. [12] 
the relative fluctuations Auq/uq = 0.05 reduce the angu- 
lar width of the transmission spectrum more than four 
times (see also Fig.[T3|). That high sensitivity makes such 
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FIG. 13: (color online) Half-width Af3 of the angular spec- 
trum for a (iii) structure. Red thin line corresponds to the 
half- width A(3 given by Eq. (ITqI) . 

a system a good candidate for use in electronic circuits 
capable of tuning the direction of charge flow. Another 
striking property is that after this abrupt drop in the 
transmission, it (i.e., the transmission) becomes practi- 
cally independent of the strength of disorder in a rela- 
tively large range, as shown in Fig. [131 

The propagation of light in analogous L-R and R-R 
disordered dielectric structures demonstrates completely 
different behavior. As the degree of disorder (variations 
Auj of the refractive indexes Uj) grows, the averaged 
angular spectra quickly reach their asymptotic "rect- 
angular" shape: a constant transmission in the region 
where all interfaces between layers are transparent fol- 
lowed by an abrupt decrease in transmission in the re- 
gion of P where the total internal reflection appears (see 
Fig. [Hbx). The frequency dependences of the transmis- 
sion coefficient and localization length have been studied 
in Ref. m. 
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FIG. 14: Mean transmission coefficient (T(/3)) for: (a) disor- 
dered type (ii) graphene structure; (b) disordered L-R dielec- 
tric structure; (c) disordered R-R dielectric structure. Differ- 
ent curves correspond to the different values of the strength 
of the disorder Auq/Uq in (a), 6n/no in (b) and (c), which 
increases from right to left from 0.1 to 1.0 with a step of 0.1. 



V. ANALYTICAL STUDY 

The features presented above can be explained, both 
qualitatively and quantitatively, in the framework of a 
rather simple theoretical approach. It can be shown^i 
that in the short- wavelength limit, kS ^ 1^ the mean 
amplitude transmission coefficient, (T*^^^), of a sample 
built of N layers, is approximately equal to 



(14) 



where tjj-^i are statistically independent complex trans- 
mission coefficients of the boundaries between the j-th 
and {j + l)-th layers, and 



It.- 



2 2cos^6>j+i 



(15) 



^■•^■+" l + cos(^,- -^,+i)' 
At small ^ <C 1, equation ([T5|) becomes 



\t,j+^f 1 - ^e]^, + -fi - (16) 



9 



arctan 



/?2 



and from Eqs. (p!4|) and (p!6|) it follows that 



1 



1 



1 



N 



(17) 



In an initially periodic array of alternating p-n and n-p 
junctions, Uj = — i^j+i = /7o, at e = [structure (iii)], 
Eq. (dll) yields 



where 



1 



1 /32 



2 (A/?)2 



In 2 



1/2 



(18) 



(19) 



is the half-width of the angular spectrum, defined as 
the value of (3 where {T^^\(3)) = 1/2. Equations (fTSj) . 
(p!9|) fit well the numerical results presented in Figs. [121 
[T3l In particular, they describe the numerically observed 
quadratic dependence on P and surprisingly weak depen- 
dence of the mean transmission on the strength of the 
disorder. 

In a disordered graphene superlattice consisting only 
of n-n and p-p junctions [structure (ii)], the mean trans- 
mission coefficient ^T^^^) at {8v?) /U^ <C 1 is given by 



N 



^0 



N 



(20) 

In this case, it is easy to see that the half-width of the 
angular spectrum strongly depends on the strength of the 
fluctuations and decreases with increasing {5u^)/Uq^ as 
can be seen in Fig. [T4k and in the relation: 



(,5w)2 



1/2 



(21) 



Note, that in the cases (i) and (iii), the transmittance 
spectrum has parabolic shape for small angles of inci- 
dence, {T^^\f3)) ^ 1 - (3^/A(3^, and the spectrum half- 
width A/3 decreases as 1/\/N when the number N of 
layers (the sample length L) increases. The same spec- 
trum property for case (ii) has been predicted in Ref. ij. 

In contrast to the charge transport in disordered 
graphene superlattices described above, the propagation 
of light in randomly layered dielectrics is similar (at 
<C 1) for L-R and RR arrays of layers with equal 
impedances (there are the analogs of p-n and p-p junc- 
tions, respectively). This follows from the fact that in 
both cases the small-angle asymptotics of the mean trans- 
mission coefficient through a boundary between layers are 
identical and at ^ <C 1 have an universal form (compare 
with Eqs. [171 [ig): 



92 



which yeilds 



(22) 



As in periodic systems, the difference in the transmis- 
sion spectra of disordered graphene and dielectric sam- 
ples [compare Eqs. (fT4j) . ([20|), and ([22|)] is a consequence 
of the above mentioned absence of imaginary part in the 
transfer matrix JVl , Eq. (|9|) . Examples of the numerically 
calculated (with no approximations) angular spectra of 
the transmission of light are shown in Fig. [T4bx. 



VI. GEOMETRICAL DISORDER 

In Section IV, we studied spatially-periodic layered 
graphene structures, in which the values of the applied 
potential in each layer were statistically independent ran- 
dom numbers. Further numerical calculations show that 
the main features of the transport and localization of 
charge in disordered graphene superlatices are rather uni- 
versal, i.e., independent of the type of disorder. When 
instead of the amplitude of the potential, the size of each 
layer is fluctuating, 

/^Uj =0, 5j = 5^ Aj 

all results are similar, at least qualitatively. In Fig. [151 
the angular dependences of the mean transmission coef- 
ficient, (T(/3)), are plotted for different strengths of the 
geometrical disorder (different values of A/S) in the case 
when Uj < e and assuming that Aj are independent 
random and homogeneously distributed in the interval 
(—A, A). As it is in the corresponding case (i) of Section 
IV (Fig. [9]), small disorder destroys the band structure of 
the underlaying periodic system, and with A/S increas- 
ing {T{/3)) takes, in the vicinity of the normal incidence 
(small /3), a parabolic form, which remains unchanged 
when increasing disorder. When e = and Uj is a peri- 
odic set of numbers with alternating signs (an array of 
random p-n junctions corresponding to the case (iii) of 
Section IV), the shape of {T{/3)) is also parabolic with 
the half- width similar to that shown in Fig. [131 



VII. CONCLUSIONS 

We have studied the transport and localization of 
charge carriers in graphene superlattices produced by ap- 
plying periodic and disordered potentials that depend 
on one coordinate. Simultaneously, the optical proper- 
ties of analogue dielectric structures composed of tra- 
ditional (right-handed, RH) dielectric and left-handed 
(LH) metamaterial layers was considered and compared 
with the charge transport in graphene. It was shown 
that in the Kronig-Penny-type periodic structures, a sort 
of total internal reflection can occur. In the case of a 
non-symmetric periodic array of alternating p-n and n- 
p junctions, along with the conduction bands and band 
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FIG. 15: (color online) Mean transmission coefficient (T(/3)) 
for graphene with geometrical disorder (section VI). Arrows 
with numbers mark the strength of disorder, A/ 6. 



Away from the normal incidence, ID random graphene 
systems manifest all features of disorder-induced strong 
localization. There exist a discrete random set of angles 
(or a discrete random set of energies for each given angle) 
for which the corresponding wave functions are exponen- 
tially localized. Depending on the type of the unper- 
turbed system, the disorder could either suppress or en- 
hance the transmission. The transmission of a graphene 
system built of alternating p-n and n-p junctions has 
anomalously narrow angular spectrum; is extremely sen- 
sitive to fluctuations of the applied potential; and, in 
some range of directions, it is practically independent of 
the amplitude of fluctuations of the potential. Our nu- 
merical results fit well the analytically calculated short 
wavelength asymptotics of the mean values of the cor- 
responding transfer matrices. The main features of the 
charge transport in graphene subject to a disordered po- 
tential have been compared with those of the propagation 
of light in inhomogeneous dielectric media. This compar- 
ison has enabled better understanding of both physical 
processes. 



gaps in the angular domain, there are also a discrete set of 
directions, in which the structures are resonantly trans- 
parent. In symmetric {ui = —1^2, £ = 0, and ^1 = ^2) 
systems, the conduction zones disappear, and the angu- 
lar spectrum of the transmission coefficient represents a 
discrete set of resonances, similar to the resonances in the 
symmetric (ni = —712, ^1 = (^2, and Zi 7^ Z2) periodic 
alternating RH-LH dielectric structures. These features 
make the direction of the charge flux easily tunable by 
changing the applied voltage. The numerical experiments 
have shown that relatively weak disorder can drastically 
change the transmission properties of the underlaying pe- 
riodic configurations. In the direction orthogonal to the 
layers created by a ID random potential, the eigenstates 
are extended for all energies and a minimal conductivity 
remains non-zero, no matter how strong the disorder is. 
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